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O . Abstract 

The physically correct parametrization of the energy levels of transition ions in crystals in terms of 
, crystal-field (CF) Hamiltonians Hqf — J2k X]g B^qC'^^ has to be based on the CF parameters Bkq 

. that lead to the correct CF splitting second moments, both the global one a and the partial ones Cfc. 

CO ' Only such parametrizations correspond to the appropriate multipolar structure of the surrounding 



CF. Each parametrization being characterized by its own multipolar crystal- held strengths Sk 
/ \ 1/2 



2k+ 



qrr l^fcgP I I foi' fc = 2,4 and 6, yields a definite second moment cr, which can be derived 
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from the additivity relationship ct^ = X^fe "'fc and the known asphericities (5'||C''^'||\I') of the central- 
ion eigenfunctions ^. The condition a = cToxp must be satisfied to ensure the parametrization's 



o 

> 

• 1—1 

X. 

^ , correctness. However, our survey of literature indicates that there exists many other well-fitted 

crystal-field parameter sets that do not obey this condition. Therefore, such sets are erratic and 
non-physical, and should be re-examined or rejected. Having a for several (> 3) eigenstates \^) 
along with the relevant (4'||C'-'^-'||^') asphericities, one can estimate Uk and Sk, which are well- 
founded experimentally. The above findings set up the parametrization process properly. Lack of 
consistency between the second moments representing various parametrizations and the pertinent 
second moments observed in experiments is presumably the main reason for deluge of formally 
accurate but accidental and inequivalent parametrizations. 

PACS: 71.15.-m, 71.23.An, 71.70.Ch, 75.10.Dg 
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1. Introduction 



The review of the crystal-field (CF) Hamiltonian HcF = J2k J2q BkqCq parametrizations available 
in literature provides a vast amount of discouraging examples of their inconsistency and divergence 
[1-6]. Only a slight part of them follows the standardization procedures [2-7]. Also the number of 
the inconsistent CF parametrizations seems to be alarming in the light of actual multipolar structure 
of a typical crystal surroundings. A suspicion arises that the reason for such widespread divergences 
among parametrizations is more fundamental than the choice of the reference system. 

Although a complete mutual orthogonality of the component TicF multipoles and, in consequence, 
their independent and separate contribution to the global second moment a have a crucial physical 
meaning, this feature is not reflected in the conventional fitting procedures. The multipolar indepen- 
dence manifests itself most fundamentally in the additivity of the relevant second moments = J2k ^k 
for any free-ion state CF splitting, where ak is a partial (multipolar) second moment produced by the 
2'^-pole T-LcF component [8-11] (section 2). Both a and corresponding to a given parametrization 
can be directly calculated. On the other hand, a and ak can be derived also from experimental data 
(section 4). There is a question - how do the two sets of values compare? 

To answer this question we have analyzed thoroughly four splitting diagrams for the states ^19/21 
^Iii/2,^Ii3/2 and ^Ii5/2 of the term of Nd^"*" ion doped into single crystal of cubic yttrium oxide 
(Y2O3) [12] (section 4). Having obtained the second moments from the splitting diagrams as well 
as the asphericities Ak{\^)) = (^||C*-'^')||^) of the considered eigenstatcs ["f) [13], and employing the 
partial second moments additivity, we have estimated the physically well-founded ak values and, in 
consequence, the relevant Sk- Unexpectedly we have found that the proposed parametrization [12] 
leads to the second moments that are incompatible with (Jexp derived from the experiment. Therefore, 
the parametrization [12] has turned out to be an accidental and non-physical one. It gives reasonable 
grounds for believing that there may exist different and inequivalent, although formally accurate, HcF 
parametrizations leading to different second moments. The erratic character of such fittings stems 
from a great number of well-fitted crystal-field parameter (CFP) sets admissible when the relevant 
second moments are not correctly fixed. 

This finding has prompted us to formulate the following criterion (the necessary condition) for the 
physical fidelity of Hcf parametrizations: each and every correct CF parametrization has to 
yield the second moment of any relevant splitting equal to the observed one. Disregard 
for this condition is presumably the main reason for an excessive number as well as inexplicable dis- 
crepancies in parametrizations available in the literature. The point is that only the fitting procedures 
which lead to cr = cJexp can be treated as potentially correct ones. 

The CFPs determination is inevitably connected with the accuracy of the initial approximation 
of the involved eigenstatcs (section 3). The eigcnfunctions applied in the presented paper have been 
obtained by simultaneous diagonalization of the interaction matrix comprising Coulomb repulsion and 
spin-orbit coupling using M. Reid f-shell programs [14]. Such type of functions (Table 2) are sufficiently 
accurate for most of the states from lower parts of the rare earth ion spectra. 

The effective character of CFPs in relation to the assumed basis of the initial eigcnfunctions 
forming the interaction matrix is shortly outlined in section 3. Possible effects of extending the applied 
approximation of the initial eigenstatcs for highly excited states is also discussed. Since these states 
might be strongly mixed by crystal-field [12,15, 16] the calculated asphericities Ak can be markedly 
changed. 

2. Basic formalism 

(k) 

The effective CF Hamiltonian HcF = J2k Y^q BkqCq for electrons with the quantum number / < 3 
is limited to the three or less 2'^-poles for = 2,4 and 6. These three terms are entirely independent 
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from each other as they transform themselves according to different irreducible representations V^'^^ 
of the three -dimensional rotation group R3. 

Each individual 2'^-pole contributes separately to the total splitting of any state. Undoubtedly, 
the relationship between a partial (rank k) second moment a^, generated by a given 2'^-polc, and 
the global second moment a, produced by the sum of the three multipoles, constitutes the adequate 
measure of the multipolar contribution [8-11]: 

-V) = h = ^E^i = j:ai (1) 

n k k 

where the gravity center of the CF sublevels within the state | J) is given by J)) = 2J+1 S with 

n 

( 

En as the energy of \n) sublevel, Sk = I 2feq:T S l-^fc^P 1 stands for the conventional CF strength 

of the 2'^-pole = Y^qBkqC^q^ [8-10,13,17], (J||CW||J) = Ak{\J)) is the dimensionless scalar 
dependent only on the angular distribution of the jJ) state electron density and reflects its 2'^-pole 
type asphericity [13], whereas \J) denotes in a brief way any free-ion initial eigenstate in approximation 
that preserves its quantum number J. Generally, when J ceases to be a good quantum number the 
state I J) is replaced by |*) (section 1). 

The CF strength Sk is an invariant of the R3 group and has clear physical meaning due to its 
direct connection with fx^ (Eq.(l)). The square of the global second moment o"^ is a simple sum of 
(t|, and any compensation of the partial contributions is excluded. The parameter S = {Yl,k ^'k)^^'^ 

or 5* = (3!^/;. 5'^) described earlier in the literature [8,17] as the global CF strength is another 
invariant of the R3 group. However, its physical meaning is rather obscure due to the scalar product 
TkSlAlin Eq.(l) [13]. 

Although the additivity of the second moments ak is crucial, this property seems to be underes- 
timated or even ignored in literure. The square of the global second moment of any CF splitting 
can be always written in the form of a positive definite sum: 

'''^1'^^) = 2JTT + SlAlm + SlAlm) . (2) 

Hence, knowing the energy spectra of n states of a paramagnetic ion embedded in a given crystal- 
field, we have at our disposal a set of n linear equations as Eq.(2). Their algebraic nature admits 
only positive solutions what markedly confines their solvability. In turn, based on the asphericities 
Ak we can determine all three Sk or at least a relation between them, depending on the number of 
available equations. All of these Sk data well-founded experimentally should be built into the applied 
parametrization procedure. 

The above linear equations (Eq.(2)) become simplified when some Sk or Ak vanish. Vanishing of 
Sk is governed by the symmetry and geometry of the central- ion surroundings. All crystallographic 
point-symmetry groups, except those of the cubic system, admit the three 2'^-poles for A; = 2, 4 and 6. 
Only in the cubic system the quadrupole term is excluded, i.e. 52 = 0. Nevertheless, Sk may vanish 
accidently due to compensation of the CFPs what results from a particular surroundings geometry. 
In turn, vanishing of some Ak depends on J quantum number representing the resultant angular 
momentum of the central ion eigenstate according to the triangle conditions for the three numbers 
(J, J,k). If J > 3, any | J) eigenstate is nominally characterized by three non-zero Ak parameters for 
k = 2,A, 6. Nevertheless, it may happen that due to the mixing of the states some Ak can compensate 
(or enhance) themselves to a different degree (Table 2). For J = 2 or 5/2 only A2 and A4 differ from 
zero leading to Eqs (2) with at most two unknowns S2 and S^- For J = 1 or 3/2 solely the quadrupole 
term is effective, i.e. only A2 ^ 0. The total splitting AE of the states in the last case allows us to 
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estimate 5*2 directly knowing accurately the relevant ^2- As an example, for any state | J = 3/2) we 
obtain = j (AEf = jS^Aj, and hence 

AE 

Since AE depends exclusively on S2, the pertinent individual CFPs cannot be determined in this 
way. 



3. The effective character of crystal-field parameters 

The actual CFPs should characterize only the multipolar structure of the central-ion surroundings in 
crystals. The point is, however, that they are always accessible in products with the corresponding 
multipolar components IC*^*^^ [|^) of electron density distribution in the central-ion open shell. Thus, 
CFPs that precisely reproduce the actual energies quite naturally depend on the ^ eigenfunction basis. 
This is why CFPs may take different values according to the basis of the assumed eigenfunctions Let 
us emphasize that only the accurate eigenfunctions provide the right CFP set which has clear physical 
meaning and a universal character as it is irrespective of a central-ion. In fact, one should realize that 
we always deal with more or less "dressed" CFPs, whereas our aim are the "bare" ones. Obviously, 
even the best fitted CFP set obtained for not particularly accurate central-ion eigenfunctions cannot 
be recognized as a definitive solution. 

CFPs obtained in fitting procedures are only roughly approximated for Russell- Saunders (RS) 
eigenfunctions, better after simultaneously taking into account Coulomb repulsion and spin-orbit cou- 
pling, and the best approximated after simultaneous diagonalization of Coulomb repulsion, spin-orbit 
coupling, and CF interactions. Moreover, within a given electron configuration, contributions of the 
above mentioned interactions to the cigcnstatcs of different energy may differ substantially. 

In the energy spectra of tripositive rare earth ions the low lying states are quite well approximated 
considering simultaneously Coulomb repulsion and spin-orbit coupling. They are somewhat worse 
approximated including the spin-orbit coupling in a perturbational way (J-mixing approach) as well 
as even for the pure RS coupling. Nevertheless, for higher excited states all these approximations 
fail. Due to the higher density of these states, i.e. proximity of their energies, the CF mixing causes 
noticeable or even considerable changes in their compositions and energies [12,15]. 

Applying uncritically some of the above simplified approximations for the whole basis of 4/^ states 
may hinder appropriate parametrization. As a matter of fact, this approach might be adequate for 
a part of lower levels, whereas for the higher ones not necessarily. The incompatibility between S'2 
estimated from the CF splitting of the '^\F^/2) state of Nd^"'":Y203 [12] (Eq.(3)) and its value resulting 
from the second moment analysis (section 4), may serve as an example of such disparity. The splitting 
AE found experimentally amounts to 196 cm~^. Using the approximate A2 value for the state ^[^^3/2) 
equal to 0.3561 (Table 2) one gets S2 = 550 cmr^, which reaches only 56% of the needed magnitude 
(section 4). It results from an overestimated A2 value, which in fact should be close to 0.200. 

Another example is an unexpected dominance of the CF splitting of |^F2) state (at 5200 cm~^) over 
that of |3P2) state (at 22750 cm ^) of Pr^"*" ion LaFs matrix [15]. It reflects well the scale of changes 
in their eigenfunctions and asphericities (mainly for the |^P2) state). For the |^F2) state we have 
AE = 143 cm~^ and a = 55.1 cm~^, whereas for the |^P2) state AE = 128 cm~^ and a = 45.2 cm~^ 
[15]. However, \Ak\ of the state |^F2) calculated within the standard approximation are less than those 
for the PP2) state and amount to: IA2I = 0.3400, |^4| = 0.0909 versus |^2| = 0.8563, |^4| = 0.1112. 

Let us remind also that the conventional CFPs contain intentionally the factor (r'^). This quantities 
- the mean k powers of the radii of the magnetic electrons [18] - concerns the central-ion itself. 
Nevertheless, since the (r^) are only simple scaling factors, it is easy to renormalize the conventional 
CFPs to the ones independent of a central-ion. 
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4. Results 



The method of adjustment of the second moments presented below has been based on the CF sphtting 
diagrams of the four lowest states originating from the term of 4/^ configuration of Nd"^+ ion in 
Nd'^"'":Y203 [12]. The Nd'^"'" ions occupy the C2 lattice sites and all the free-ion states are split into 
Kramers doublets. The data come from the optical absorption and fluorescence spectra reported by 
Chang et al [12] (Table 1). The relevant asphericities Aj. of the involved states calculated for the 
eigenfunctions obtained using M.Reid f-shell programs [14] and the free-ion data reported by Carnall 
et al [16] are compiled in Table 2. 

The additivity of the partial second moments leads to the following set of linear equations for 
(Table 1): 

1 

10 
1 
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1 

14 
1 
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SlA\ 




]+SlAl 


("^^3/2) + S^Al 


(^^3/2) 


SlAl 


('/l5/2^ 


)+SlAl 


{%./2)+SlAl 


^15/2): 



^15/2) . 



(4) 



The RS labels in the parentheses arc assigned according to the maximal component in the free-ion 
eigenvector of the state. Substituting the values of i^Ij) from Table 1 and A^ (^Ij) from Table 2 
into Eqs (4) we obtain in [(cm~^)^]: 



^ (0.22645| + 0.224951 + 1.16065'^ 
^ (0.2478S| + 0.15185| + 0.1138S|) 
^ (0.31535| + 0.22095| + 0.38275^) 
(0.4339S| + 0.47425| + 3.1560S|) 



60238 
33507 
40373 
83211. 



(5) 



Si = 



950000 {cm-^f 



Si = 



We have found some acceptable solutions centered at the values 
950000 (cm-i)2 and = 150000 {cm-'^f which can be recognized as reliable. This non-optimized 
solution satisfies Eqs (5) with rms deviation 368 (cm~^)^, i.e. within the accuracy of ca. 99%. 

On the other hand the CFPs for Nd2+:Y203 given in Table VI [12] yield the following in 
[{cm-^f]: Si = 278472, Sl = 756148, and S^ = 76576, which are smaller than those calculated from 
Eqs (5), particularly Sl- Substituting now the asphericities Ak of the involved states (Table 2) into 
Eq.(l) we obtain the following cj^ in [(cm~^)^]: 33094, 16548, 20926 and 46250 for the corresponding 
split states ^-^11/2) ^-^13/2 aiid ^/i5/2- They reach only a half of the observed (Table 1). So, the 

proposed parametrization does not guarantee the required second moments of the observed splittings, 
and in consequence is devoid of physical meaning. Its pretty high accuracy {rms = 7.0 cm~^) is 
misleading as it refers to an accidental and improper a. The correct parametrization has to provide 
any other CFP set that should ensure both the same energy of the CF sublevels and the correct second 
moments. Such solution is always provided by each equation (t^(| J)) = J2k SlAk'^{\J)). 



5. Discussion 

A typical example of CF parametrization analyzed in section 4 reveals that the set of CFPs derived 
from spectroscopic data in an advanced fitting procedure (Table VI in [12]) yields only about 50% 
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of the observed squares of the global and partial second moments of the splittings. It means that 
the physically correct parametrization, which gives the appropriate second moments, turns out to be 
correlated with a quite different a = (Texp- It points to the important conclusion that many coincidental 
T-LcF parametrizations can be found for any considered splitting if there are no restrictions imposed 
on the relevant second moments. Such possible accidental parametrizations can differ in the following 
ways: 

(i) in the global second moment ct, 

(ii) in a different partition = J2k same a, 

(iii) in a different mutual spatial orientation of the component multipoles [19] at the same a and the 
same ak [19]. 

Hence, all the CFP sets fulfilling conditions (i)-(iii) are, to a different extent, incquivalent and therefore 
false. The correct and equivalent parametrizations should yield a = (Tcxp, the identical proper partition 
£7^ = J2k "^fc the same mutual spatial orientation of the component multipoles. All the remaining 
parametrizations not necessarily correspond to certain local rms minima occurring in the fittings, but 
can be merely accidental and non-physical ones which should be verified. Due to their diversified 
Sk (the rotational invariants) or various spatial orientations of the component multipoles they do 
not follow the standardization procedures based on the reference frame rotations [2,3,4,7,19]. While 
verifying inequivalent parametrizations (concerning the same paramagnetic ion in a lattice system) 
one should a priori impose the known Uexp- It should rcdTicc the parametrizations to the correct and 
equivalent ones, differing at most in the relevant reference frame. 

Disregard for the second moments in TicF parametrization explains well the inflation of various 
CFP sets obtained in fittings. There are no other convincing reasons for such deluge of solutions. After 
all, the actual multipolar structure of the parametrized CF unknown during the fitting process is rather 
strictly defined. A question arises - what is the probability of a fictitious partition cj^ = af. in the 
fitting process when a = cxexp is fixed a priori? Such probability should not be high since each partition 
corresponds to characteristic proportions of the component multipoles in ^CF- So, a coincidence of the 
inequivalent parametrizations seems to be an exceptional case. The orientation of the nominal reference 
frames (unknown by definition) [4,6] during searching for the proper partition is inessential. Let us 
notice that identities of both a and ak correspondingly for different coincidental parametrizations are 
the necessary conditions for their equivalence and fidelity. However, some equivalent CFP sets can be 
false if their identical a and ak are not equal to the respective actual (observed) values. 

In order to derive partial a^ and hence Sk from the relevant linear equations (Eqs (4)) and known 
c"exp, we still require the asphericities of the involved states. The accuracy depends on how ade- 
quate are the assumed eigenfunctions and hence determines the precision of both a^ and S^. Therefore, 
it is recommended to use the CF splitting data for low-energy states of well-defined eigenfunctions. 
On the other hand, in some cases even the simple RS functions can turn out to be sufficiently good 
(Table 2). For example, the eigenfunctions of the states \'^H^), \^F2) and \^Pi) for configurations 4/^ or 
4/^^ in Pr^+ or Tm^+ ions remain single-component eigenvectors (within the assumed approximation) 
since in these configurations there are no other states with the same J. 

For higher excited states, to approximate the eigenfunctions sufficiently well, the full simultaneous 
diagonalization of the complete interaction matrix including the CF Hamiltonian is needed due to the 
importance of CF effect. But in such case a more sophisticated self-consistent numerical approach to 
the problem should be applied. 

The CF analysis carried out in [12] was based on a Hamiltonian of C2 point group symmetry 
including J-mixing effects. The program by M. Reid [14] employs a simultaneous diagonalization of 
the interaction matrix including both Coulomb repulsion and spin-orbit coupling. The discrepancies 
in the resulting A^'s are negligible. Actually they are close to the asphericities of the pure RS states 
(Table 2). 
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The hnear equations (see Eq. (2)) relating and Aj^ can be employed alternatively. One can find 
Sk using Ak, or vice versa, obtain A^ from Sk- Having the global second moments a for several states, 
e.g. Eqs (4), we may reconcile these quantities as it has been shown in section 3 for the 1^-^3/2) and 
four lowest states of the Nd^+:Y203 [12]. 

6. Conclusions 

Based on Nd^^:Y203 spectroscopic data and the related CF analysis [12] we have shown that formally 
correct CF parametrization does not lead to the experimental second moments of the relevant CF 
splittings. There are reasonable grounds for believing that it is not an exceptional case among available 
CF parametrizations. Presumably, it can be also the main reason for their inconsistencies. Thus, 
the fundamental criterion to be fulfilled by the physically founded and correct parametrization has to 
yield the second moment of any pertinent splitting equal to the observed one. This necessary condition 
sounds trivially only because there exists a widespread believe that it is satisfied automatically in used 
fitting procedures. In order to verify a given parametrization it is enough to know the asphericities 
= (^'||C(^)||^'), = 2,4 and 6, for the initial central -ion eigenfunctions ^' along with the 
corresponding Uexp and check the relationship a^^^ = J^k '^k^ti'^)- avoid false, non-physical 
parametrizations this condition should be explicitly built in the fitting programs. 
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Table 1: Experimental crystal-field splitting data for states of Nd^''":Y203 [12] 



State * 


Centroid parameter 


Number of 


Total splitting 


Square of the second 


[SL]J 


[cm~^] 


CF doublets 


AE [cm-i] 


moment o"^ [cm~^]^ 


h/2 


314 


5 


643 


60238 


hl/2 


2171 


6 


462 


33507 


4t 

■■-13/2 


4133 


7 


515 


40373 


%5/2 


6160 


8 


770 


83211 



* Russell- Saunders labels according to the maximal component in the free- ion eigenvector of the state 



9 



Table 2: Asphericities of chosen Nd^"*" free-ion eigenstates (the values for the upper component of the 
state are given in parentheses) 



No 


*Eigenstate \^) 
Energy [cm~^] 


Multipolar asphericity 

A2{^) A4(*) Aq{^) 

(*||C(2)||^-) (*||C(^)||*) (*||C(6)||*) 


1 


0.9844|4l9/2)-0.1647|^H(2)9/2)+0.0566|^H(l)9/2) 

+ 2 components 




-0.4758 


-0.4742 


-1.0773 




(-0.4954) 


(-0.4904) 


(-1.1085) 


2 


0.9947|nn/2)-0.0949|^H(2)n/2)+0.0361|^H(l)n/2) 

+ 1 component 

1862 

I^Ill/2) 


-0.4978 
(-0.5045) 


-0.3896 
(-0.3935) 


-0.3374 
(-0.3399) 


3 


0.9979|ni3/2)+0.0609|^Ki3/2)-0.0231|^Ii3/2) 


-0.5615 


-0.4700 


-0.6186 




3845 


(-0.55G9) 


(-0.1(391) 


(-0.(3217) 


4 


0.9938|4li5/2)+0.1110|^Ki5/2) 


-0.6587 


-0.6886 


-1.7765 




5907 

I^Il5/2) 


(-0.6438) 


(-0.6850) 


(-1.7999) 


5 


0.9698|4F3/2)+0.2232|^D(l)3/2)-0.0614|^D(2)3/2) 

+ 3 components 

11381 

l%/2) 


0.3561 
(0.3578) 




(0) 




(0) 


6 


0.9879|4F5/2)+0.1468|^D(l)5/2)-0.0321|^F(2)5/2) 

+ 3 components 

12420 

"^^5/2) 


0.3277 
(0.3220) 


0.1843 

(0.1890) 




(0) 


7 


0.9648|*F7/2)+0.2004|^G(l)7/2)-0.1582|^G(2)7/2) 

+ 3 components 

13383 


0.4650 
(0.4601) 


0.0076 
(0.0537) 


-0.4233 
(-0.4552) 


8 


-0.8670|4F9/2)-0.4420|^H(2)9/2)+0.1526|^H(l)9/2) 

+ 4 components 

14652 

l%/2) 


0.5748 
(0.7136) 


-0.3050 

(-0.4051) 


0.1040 

(0.1799) 


9 


0.7205|^P3/2)+0.6371|^D(l)3/2)-0.1869|^D(2)3/2) 

+ 3 components 

26179 

PP3/2) 


0.2818 
(0.2981) 




(0) 




(0) 



*The eigenfunctions and eigenvalues have been calculating using M. Reid f-shell programs [14] and free-ion data 
reported by Carnall et al [16]. Only the first three main components of the states are given explicitly. 
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